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ABSTRACT 

We investigate the gas content and baryonic Tully-Fisher relationship for 
extremely low luminosity dwarf galaxies in the absolute magnitude range —13.5 > 
M r — 51og 10 /i7o > —16. The sample is selected from the Sloan Digital Sky 
Survey and consists of 101 galaxies for which we have obtained follow-up HI 
observations using the Arecibo Observatory and Green Bank Telescope. This 
represents the largest homogeneous sample of dwarf galaxies at low luminosities 
with well-measured HI and optical properties. The sample spans a range of 
environments, from dense groups to truly isolated galaxies. The average neutral 
gas fraction is (/ gas ) = 0.6, significantly exceeding that of typical gas-rich galaxies 
at higher luminosities. Dwarf galaxies are therefore less efficient at turning gas 
into stars over their lifetimes. The strong environmental dependence of the gas 
fraction distribution demonstrates that while internal processes can reduce the 
gas fractions to roughly / gas = 0.4, external processes are required to fully remove 
gas from a dwarf galaxy. The average rotational velocity of our sample is (v Iot ) = 
(W20j ) t/2) = 50 km s -1 based on HI line-widths. In this luminosity range, the 
optical Tully-Fisher relationship has significantly more scatter compared to the 
baryonic relationship. Including more massive galaxies from the literature, we fit 
a baryonic Tully-Fisher slope of M baryon oc v^ t 0±0A5 . This slope compares well 
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with CDM models that assume an equal baryon to dark matter ratio at all masses. 
While gas stripping or other processes may modify the baryon to dark matter 
ratio for dwarfs in the densest environments, the majority of dwarf galaxies in 
our sample have not preferentially lost significant baryonic mass relative to more 
massive galaxies. 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics 

1. Introduction 

The properties of low mass dwarf galaxies, and their differences from luminous galaxies, 
provide important clues to understanding both cosmology and galaxy formation. As probes 
of galaxy formation, dwarf galaxies should be extremely sensitive to processes such as star 
formation feedback, ram pressure stripping, and tidal stripping. These processes are often 
invoked to explain properties of relatively luminous galaxies and should even more dramat- 
ically affect dwarf galaxies, which have smaller masses and shallower potential wells. As 
probes of cosmology, the number of dwarf galaxies has provided a challenge for cosmological 
models. 

The current Cold Dark Matter model with a cosmological constant (ACDM; Spergel 
et al. 2006) makes robust predictions for the number of dark matter halos as a function 
of mass (Sheth et al. 2001). However, extending this model to predict the number density 
of observed galaxies requires including physical effects such as gas cooling, star-formation, 
supernova feedback, etc., that are too complex to follow exactly in numerical predictions. It 
is clear that the number of galaxies as a function of luminosity is much shallower than the 
prediction for the number as a function of halo mass (Benson et al. 2003). The most extreme 
example of this problem is the "substructure problem" in the Local Group (Klypin et al. 
1999; Moore et al. 1999). If we assume galaxies are associated one-to-one with dark matter 
halos and that ACDM is correct, the mass-to-light ratios of galaxies therefore must decline 
as a function of host halo mass. That is, the efficiency of forming stars and/or retaining gas 
must be much smaller for low mass galaxies. 

A number of physical models have been proposed to explain high mass-to-light ratios 
in low mass galaxies, including photoevaporation of gas during reionization (Shapiro et al. 
2004; Mayer & Moore 2004), suppression of gas accretion during reionization (Gnedin 2000; 
Bullock et al. 2000; Benson et al. 2003), feedback from supernovae (Dekel & Woo 2003; 
Nagashima & Yoshii 2004), and ram pressure stripping of gas (Mori & Burkert 2000). Some 
of these processes only occur for dwarf galaxies near a luminous galaxy and not for isolated 
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ones (Kravtsov et al. 2004). The relative importance of these processes in the evolution 
of dwarf galaxies has significant implications for galaxies of all sizes at both low and high 
redshift. For example, the importance of galactic winds in low mass galaxies will strongly 
affect our understanding of the mass-metallicity relation (Tremonti et al. 2004; Lee et al. 
2006) and the enrichment of the intergalactic medium at high redshift (Scannapieco et al. 
2005). By comparing the gas content and Tully-Fisher relationships of dwarfs as a function of 
environment, it should be possible to determine the importance of "external" environmental 
effects like stripping relative to "internal" effects such as feedback. 

Historically, extremely low luminosity nearby galaxies have been difficult to find because 
their colors and apparent sizes are similar to those of more luminous, and far more numerous, 
background galaxies. By far the best known set of dwarfs is that in the Local Group (Mateo 
1998; Grebel et al. 2003; Dolphin et al. 2005), where the even the lowest mass dwarfs can 
be detected as over- densities in stellar counts (Willman et al. 2005; Zucker et al. 2004). The 
dwarf galaxy population has also been well studied in nearby clusters, such as Virgo and 
Fornax, where the number of dwarfs is large relative the background population (Binggeli 
et al. 1985; Ferguson & Sandage 1991) or in other environments where dwarfs are expected, 
eg. as satellites to bright galaxies (Zaritsky et al. 1997; Gutierrez et al. 2002). Assembling 
samples of dwarfs in more distant and less dense environments has been challenging. Using 
criteria such as surface brightness and irregularity, several authors have assembled large 
catalogs of nearby dwarf galaxies (Impey et al. 1996; Schombert et al. 2001; Hunter & 
Elmegreen 2004; Karachentsev et al. 2004). 

Addressing the above issues requires a large sample of galaxies across a range of en- 
vironments, with minimal selection effects. Deep, wide-field optical spectroscopic surveys, 
such as the Sloan Digital Sky Survey (SDSS), provide a large, homogeneous sample of dwarf 
galaxies covering a wide range of galactic properties and environments. Concurrently, blind 
HI surveys such as the HI Parkes All Sky Survey (HIPASS; Meyer et al. 2004) and the up- 
coming Arecibo Legacy Fast ALFA Survey (ALFALFA; Giovanelli et al. 2005) will provide 
similarly large samples, with a very complementary set of selection effects. 

Unlike present-day massive disk galaxies, which have used up much of their gas reser- 
voirs (Read & Trentham 2005), the gaseous component of dwarf galaxies cannot be ignored 
(Schombert et al. 2001; van Zee 2001). Indeed, we show that at the lowest luminosities 
dwarf galaxies are mostly neutral gas. A galaxy's gas content reflects the efficiency of star 
formation modified by processes which add or removal gas from the system such as gas in- 
fall or ram pressure stripping. The goal of this paper is to combine the optical properties, 
dynamics and neutral gas content of low mass dwarf galaxies in order to better understand 
the processes important in dwarf galaxy evolution. 
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We present neutral hydrogen (HI) observations for 101 low luminosity dwarf galaxies 
selected with uniform and well-understood criteria from the SDSS. We describe the catalog 
from which these dwarfs were selected in §2.1, estimate the environment of each object in 
§2.2 and discuss the HI observations in §2.3. We investigate correlations between the gas 
fraction, color, star formation rate and environment of this sample in §3.1 and the Tully- 
Fisher relationship in § 3.2. When necessary, we assume a Hubble constant of H = 70/i7o Mpc 
km s -1 where hj = 1. 



2. Data 
2.1. The Sample 

We selected dwarf galaxies with magnitude M r — 51og 10 /i 70 > —16 for HI observations 
from the low luminosity galaxy catalog of Blanton et al. (2005), a subsample of the New 
York University Value-Added Galaxy Catalog (NYU-VAGC; Blanton et al. 2005b) based on 
the SDSS (York et al. 2000) Data Release 2 (DR2; Abazajian et al. 2004). The Blanton 
et al. (2005) catalog represents a significant improvement over simply selecting galaxies from 
the SDSS catalog, which is not optimized for nearby, low surface brightness galaxies. It 
is a flux-limited sample selected with an optical apparent magnitude limit of m r ~ 17.8. 
For each galaxy, the catalog provides the SDSS redshift, multi-band photometry, structural 
measurements (in particular, two-dimensional Sersic fits as described in Blanton et al.) and 
emission line measurements (Tremonti et al. 2004). Distances are estimated based on a 
model of the local velocity field (Willick et al. 1997). Distance errors have been folded into 
error estimations of all distance-dependent quantities such as absolute magnitude and HI 
mass. 

There are two selection effects of the Blanton et al. (2005) catalog relevant to this 
analysis: (1) The SDSS flux limit at m r ~ 17.8 restricts the detection of galaxies with 
M r — 5 log 10 h 70 = —16 to within 60 Mpc in distance (20 Mpc for the least luminous galaxies, 
M r — 51og 10 /i7o = —13.5). The catalog therefore does not span the full range of galaxy 
environments. While there are plenty of voids and groups in this volume, there are few 
very dense regions. For example, given the angular limits of the SDSS DR2, we see only the 
fringes of the Virgo Cluster. (2) The catalog does not include low surface brightness galaxies. 
As shown in Blanton et al. (2005), the SDSS completeness as a function of half-light surface 
brightness drops below 50% at fi 50 . r ~ 23.5 mag arcsec -2 . In the absolute magnitude range of 
interest, we would miss half of the dwarf galaxies in the Local Group (Mateo 1998) given this 
surface brightness cutoff. According to the definition of McGaugh (1996), 10% of our sample 
galaxies are considered Low Surface Brightness (LSB) with /if > 22.7. With the exception of 
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these two selection effects, the Blanton et al. spectroscopic catalog is not selected by galaxy 
property in any particular way. For example, galaxies have not been selected for or against 
irregular looking morphologies. 

We present here HI observations for 101 dwarf galaxies in the absolute magnitude range 
— 13.5 > M r — 5 log 10 h 70 > —16 selected from the Blanton et al. (2005) catalog. These data 
are part of a larger, on-going survey to obtain HI line-widths and Ha rotation curves for a 
significant fraction of galaxies in this catalog with M r — 51og l0 hj > —16. The complete 
catalog contains 1034 galaxies at these luminosities. The 101 galaxies presented here are 
representative of the full catalog. 

2.2. Environment and the Nearest Neighbor Distances 

A goal of this paper is to understand the effects of environment on the properties of 
dwarf galaxies. To quantify the environment, we calculate for each dwarf galaxy the projected 
distance to its nearest luminous neighbor. We cannot rely on the SDSS alone to determine 
this distance for several reasons. First, the angular distances between nearest neighbor 
galaxies can be large for this nearby sample — for example, searching a 1 Mpc region around 
a galaxy 30 Mpc away corresponds to 2 degrees on the sky. Many of our dwarf galaxies are on 
the SDSS Southern stripes, which are only 2.5 degrees wide. In addition, because the SDSS 
reduction software is not optimized for large extended objects, and fails to process them 
correctly, the SDSS catalog does not contain many of the bright galaxies within 30 Mpc. 
Thus, to calculate the environments of our dwarf galaxy sample, we need a supplemental 
catalog that extends beyond the SDSS area and contains the brightest galaxies. 

Both of these considerations drive us to use the The Third Reference Catalog of Galaxies 
(RC3; de Vaucouleurs et al. 1991), which is a nearly complete catalog of nearby luminous 
galaxies. To determine environments for our dwarf galaxies, we determine the distance to 
its nearest "luminous" neighbor. We define "luminous" galaxies as brighter than M r — 
5 log 10 hjo < —19 or have circular velocities greater than V c > 100 km s -1 (Haynes et al. 
1999). While we would like to include galaxies of all sizes in our environment estimate, 
the availability of a complete comparison catalog restricts our search to luminous neighbors. 
From the B and V photometry listed in RC3, we infer M r for each RC3 galaxy using the 
kcorrect product (see §3.1) . For galaxies which have the relevant entries listed, we call 
galaxies luminous if M r — 51og 10 /i7o < —19. While this absolute magnitude is set by the 
catalog limits, it roughly corresponds to the dividing line between normal and dwarf galaxies. 
For RC3 galaxies which do not have the relevant entries, but do have HI data listed, we call 
them luminous if W20 > 200 km s _1 (W20 is twice the maximum circular velocity of the HI 
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gas, see §2.3). Finally, there are some galaxies with neither HI data nor optical photometry 
listed in RC3. For this small set, we extract the "magnitude" from the NASA Extragalactic 
Database (NED), which empirically is very similar to the B band RC3 magnitude, and 
guess M r based on that magnitude. Additionally, we update the coordinates in RC3 using 
the NED coordinates for each of the catalog objects. This set of bright galaxies is not 
perfectly uniform, but is suitable for our purposes. 

To determine the nearest neighbor distance, we ask whether there is a luminous RC3 
neighbor within a projected distance of 2.5 Mpc and 300 km s _1 in redshift for each dwarf 
galaxy in our sample. We choose these limits because Blanton et al. (2006) have shown that 
environmental affects are strongest inside this scale. About 20% of our galaxies have no 
such neighbor; for these objects, we set the nearest neighbor distance to 2.5 Mpc. Using 
the SDSS Catalog Archive Server*, we have looked at wide field mosaic images around all 
of our objects and checked by eye that the RC3 catalog object is indeed the nearest bright 
neighbor. This check is not perfect due to the edge effects of SDSS, but provides a sanity 
check on our procedure. 



2.3. HI Observations 

We obtained HI radio observations with the Arecibo 305-m telescope and the Green 
Bank 100-m Telescope (GBT) in July - August 2005 and March - April 2006. At Arecibo, 
we used the L-wide receiver in 9-level sampling mode with 1024 channels and a 12.5 MHz 
bandwidth. At GBT, we used the L-band with 8192 channels and a bandwidth of 12.5 MHz. 
The Arecibo configuration yields a velocity resolution of 2.5 km s _1 , the GBT configuration 
yields a resolution of 0.3 km s" 1 . Standard ON/OFF position switching was used during 
both the Arecibo and GBT observations. The total on-source integration times were on 
average 10 minutes for the Arecibo data and 20 minutes for the GBT data. The average 
half-light radius of the dwarf galaxies in our sample is ~ 8" and they should be completely 
contained within the radio beamsize of 3' and 9' for Arecibo and GBT, respectively. We 
have searched the beam coverage for nearby contaminating galaxies and have removed the 
few cases in which a nearby galaxy is contained in the beam and is within 100 km s^ 1 of the 
targeted galaxy. Optical images and HI profiles for three representative galaxies are shown 
in Figure 1. Our observations yielded strong HI detections for 88 of 101 galaxies. The 13 
galaxies which did not yield HI detections are represented by upper limits in the analysis of 
§3.1 and are excluded from the Tully-Fisher analysis in §3.2, since we cannot estimate HI 
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line-widths for these galaxies. 

We compute the 20% HI line- width (W20) by finding the peak HI flux within 150 km s~ x 
of the optical radial velocity of each galaxy and computing the difference between the nearest 
points having 20% of the peak flux. The integrated HI flux is calculated by expanding the 
W20 values by 20 km s" 1 on a side and integrating the flux in this region. Errors bars on 
the line widths and integrated fluxes were computed using a Monte Carlo bootstrap method: 
noise was added to the stacked one-dimensional radio spectra based on the observed variance 
in the baseline and the observed quantities measured. Error bars on the line-width and 
integrated flux were calculated from the scatter in the mean quantities recovered form the 
Monte Carlo simulations. We convert the HI integrated flux into HI mass based on the 
optically-thin approximation: M m = 2.356 x 10 5 -D 2 F H i where D is the distance in Mpc and 
F H i is the integrated HI flux in Jy km s" 1 (Eqn. 9; Haynes & Giovanelli 1984). We note that 
this standard HI mass conversion is a lower limit to the true HI mass if the galaxy contains 
optically thick HI. Error bars on the HI mass include the HI flux error and the error in the 
distance. 



3. Results 

The 101 dwarf galaxies presented in this paper represent the largest homogeneous sample 
of galaxies fainter than M r — 5 log 10 /170 > —16 with well-measured HI and optical properties. 
Properties of the sample are listed in Table 1; the distribution of these properties is shown in 
Figure 2. The average absolute magnitude, corrected for Galactic extinction (Schlegel et al. 
1998), of our sample is (M r ) = —15.4 at an average distance of SO/i^Mpc. Galaxies in the 
sample tend to be blue, (g — r) = 0.3, reflecting the fact that our sample does not probe 
the densest environments where red dwarf ellipticals tend to reside. The galaxies have an 
average effective surface brightness of (/i r ,cff) = 22.4 mag arcsec -2 , and exponential surface 
brightness profiles as measured by the Sersic index ((n Se rsic) = 1-1)- The average effective 
half-light r-band radius is (r cff r ) = 8" = l.O/i^kpc. These galaxies exist in environments 
ranging from dense groups to truly isolated objects: about 15% of the sample are satellites 
of luminous host galaxies, defined as being within 200 kpc projected distance and 300 
km s _1 in velocity of a galaxies with M r — 51og 10 /i7o < —19 (see §2.2). Meanwhile, about 
65% of the galaxies are further than 500 kpc of any such host galaxy — they are very 
isolated galaxies. 

The dwarf galaxies in our sample were not selected with respect to morphology or gas 
content. The sample consists of ~ 10% gas-poor dwarf galaxies commonly known as "dwarf 
spheroidal" or "dwarf elliptical" galaxies. The majority of our sample consists of gas-rich 
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dwarfs commonly known as "dwarf irregular" galaxies. This latter term is commonly mis- 
interpreted to mean that gas-rich dwarf galaxies have disorganized morphology or turbulent 
kinematics. Historically, the term "irregular" indicates that the HII regions in a galaxy do 
not follow a regular spiral arm pattern as observed in normal massive disk galaxies (Sandage 
1975; Skillman 1996). Indeed, our sample dwarf galaxies do not show any evidence of spiral 
structure, based on the visual inspection of unsharp-masked SDSS r-band images. However, 
our dwarfs do have well-defined optical centers and regular morphology (Figure 1) as has 
been noted by previous studies of dwarf galaxies in this regime (Matthews & Gallagher 2002). 
Our HI radio profile shapes suggest that our dwarf galaxies have regular rotation fields: 18% 
of the HI profiles have distinct double-peaked shapes implying a coherent rotation field; this 
fraction increases to 30% if we limit the sample to edge-on galaxies (b/a < 0.4) only. Thus 
while the classification "dwarf irregular" is correct for the majority of our sample galaxies, 
due to its common misinterpretation, we simply refer to all the galaxies in this sample simply 
as low luminosity dwarf galaxies. 



3.1. HI Gas Fractions 



We first investigate the ratio of gas mass to total baryonic mass in our dwarf galaxy 
sample. We define the gas fraction as / gas = M gas /(M gas + M ste n ar ). The denominator of 
this expression is the total baryonic mass defined as the sum of the gas and stellar mass. 
We take M gas = l.AMm, where Mri is calculated in §2.3 and the multiplicative factor 
takes in account the presence of helium and metals. We do not include a correction term 
for molecular hydrogen (H 2 ) as the H 2 /HI ratio is much lower in dwarf galaxies compared 
to luminous spirals (Leroy et al. 2005; Taylor et al. 1998). This deficiency holds even after 
accounting for a strong metallicity dependence of the H 2 tracer molecule CO (Garnett 2002). 
The stellar mass, M ste n ar , is based on the optical SDSS i-band magnitude and (g — r) colors 
using the mass-to-light ratios determined from the galaxy evolution models of Bell et al. 
(2003) using the Kroupa et al. (1993) initial mass function (that is, we subtract 0.15 dex 
from Bell et al. Table 7 calculated for a 'diet' Salpeter initial mass function). We choose this 
combination of magnitude and color to minimize contamination from emission lines and thus 
probe the "true" stellar continuum (West 2005). We have calculated the stellar mass based 
on the broad-band fluxes using the kcorrect product ^ and based on the spectra published by 
Kauffmann et al. (2003). Both these methods produce somewhat lower stellar mass (~ 30%) 
than predicted by Bell et al. We stress that the mass-to-light ratios are very sensitive to 
the assumed initial mass function. We note that the gas fraction does not depend on the 
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distance to individual galaxies. 

The average gas fraction for our sample is (/ gas ) = 0.6, with gas fractions as high as 95%. 
Thus, for the majority of the dwarf galaxies in our sample, the baryonic mass is dominated 
by gas, rather than stars. In Figure 3 we plot gas fraction as a function of absolute r-band 
magnitude. The flux limits of our sample are clearly seen in this figure. We compare our 
results to representative samples from the literature by converting the literature magnitudes 
to r-band using the photometric transformations from kcorrect (see the web page referenced 
above) and recalculate gas fractions using the same methods described above. The most 
comparable sample of galaxies is that of LSB dwarf galaxies by Schombert, McGaugh, & 
Eder (2001). The gas fractions in this sample are comparable to our SDSS (higher surface 
brightness) sample. For more luminous spiral galaxies, the average gas fraction is much lower: 
Bell & de Jong (2000) find an average gas fraction of 30% for a sample of low-inclination 
spiral galaxies. The average gas fraction for the large field spiral sample of Haynes et al. 
(1999) is 25%, assuming an average color V — I — 1.0. A clear trend is seen in Figure 3: 
less luminous galaxies tend toward higher gas fractions. While galaxies exist in the lower 
left region of the plot (faint, gas-poor dwarf elliptical or dwarf spheroidal galaxies), they 
are rare in the upper right region (based on the HIPASS sample, West et al. in prep). In 
other words, galaxies fainter than M r — 5 log 10 h 70 ~ —16 span the full range of possible 
gas fractions, while there is an upper limit to the gas fractions of brighter galaxies which 
depends on absolute magnitude. This trend has been noted by several authors (McGaugh & 
de Blok 1997; Schombert et al. 2001; Lee et al. 2002); Figure 3 extends this to significantly 
fainter absolute magnitudes compared to previous studies. 

We interpret the gas fractions in Figure 3 as evidence that, compared to more luminous 
galaxies, dwarf galaxies are far less efficient at turning gas into stars over their lifetimes. 
This is likely related to lower on average gas surface densities in lower luminosity systems. 
Kennicutt (1998), Martin & Kennicutt (2001) and references therein have shown for luminous 
galaxies that star formation does not occur below a critical gas density (~ 5 M pc -2 ). Hunter 
et al. (1998) confirm that this limit holds for low mass systems, although the critical density is 
slightly lower (~ 3 M pc~ 2 ), due in part to the lack of spiral arms and solid-body kinematics 
in the inner regions of dwarfs. We cannot precisely measure the average gas density in this 
sample because our HI single-dish observations do not provide an estimate of the HI radius. 
However, we can estimate the HI radius by assuming it scales (roughly) with the optical 
radius. We adopt an average Hi-to-optical radius scaling calculated by Begum et al. (2006) 
of r m /r Ua = 2.5. The average gas density of our sample is 1.5M Q pc~ 2 with more than 
a factor of two scatter. Thus, the majority of galaxies in this sample are at, or below, 
the critical gas density required for star formation. If the threshold described in Kennicutt 
(1998) exists, the star formation in these dwarf galaxies implies that the local density of gas 
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in the regions of star formation significantly exceeds the global density averaged over the 
whole galaxy. 

In the top right panel of Figure 4, we show that the gas fraction is a strong function of 
g — r color in the sense that bluer galaxies in our sample are gas-rich. We assume that the 
spread in g — r color represents a spread in time-averaged star formation over the past 1 Gyr, 
as opposed to a spread in metallicity, for two reasons. First, the bluest colors in our sample 
cannot be explained by old, metal poor stellar population synthesis models. Second, the 
majority of galaxies in our sample are currently forming stars (as indicated by Ha emission 
in the SDSS spectra) which dominates their broad-band colors (Bruzual & Chariot 2003). 
Thus, in the top left of Figure 4, the gas-rich galaxies have had a higher fraction of their 
star-formation in the recent past. In contrast, the current star formation rate per unit area 
(within the last few million years) is not correlated with gas fraction (Figure 4, bottom left). 
We calculate this current star formation rate based on the w-band flux using the conversion 
from Hopkins et al. (2003). For comparison, we also calculated the star formation rates 
based on the Ha emission line flux in the SDSS spectra (detected for most of our galaxies), 
but found that the resulting rates suffered from the uncertainty of extrapolating the Ha 
flux from the 3" SDSS fiber aperture. We normalize the total star formation rate by the 
area over which star formation is occurring in order to remove any size-dependent variations 
within our sample. The median total star formation rate per unit area based on the w-band 
for our sample is SFR = 1 x 10~ 3 M yr _1 kpc _2 ; These rates are comparable to previous 
studies of gas-rich galaxies in this luminosity range (Hunter & Elmegreen 2004; van Zee 2001; 
Barazza et al. 2006). The star formation rates per unit area are shown in the bottom left of 
Figure 4 and are not correlated with gas fraction. The gas fraction is better correlated with 
the time-averaged star-formation rate over the last billion years (g — r color) than with the 
recent star- formation rate in the last few million years. This implies that the star formation 
rates we are observing now are different than those in the recent past; the star formation 
rates appear to be stochastic. 

We attempt to quantify the stochastic nature of star formation in our sample by calcu- 
lating the stellar birthrate timescale. In the left panel of Figure 5, we plot the observed stellar 
mass divided by the current star formation rate (M stenai / SFRi^). If the current rate of star 
formation has been constant over the lifetime of a galaxy, this quantity would represent the 
time needed to form its observed stellar component. For our sample, the median timescale to 
make stars is 15 Gyrs with a nearly symmetric scatter to both shorter and longer timescales. 
Observations of the color-magnitude diagrams of dwarf galaxies in the Local Group suggest 
that all dwarfs have an ancient stellar population component (Mateo 1998; Dolphin et al. 
2005). If we assume all the galaxies in our sample formed at early times (10-12 Gyr), this 
implies that roughly half the sample had lower star formation rates and half had higher rates 



- 11 - 



in the past than presently observed. This line of argument again suggests that star formation 
in dwarfs fluctuates and that the current star formation rate is not necessarily representative 
of the time-averaged rate. In the right panel of Figure 5 we also investigate the future of 
star formation by plotting the gas consumption timescale (M gas / SFR tot ). In this plot we 
have not included gas-recycling, which would increase the gas consumption timescales. The 
majority of galaxies in this sample will be able to make stars, and therefore consume gas, at 
their present rates for another Hubble time. 

The strongest observed trend in our sample is that between the gas fraction and g — r 
color (Figure 4, top left panel). This trend may be due to either internal or external processes. 
For example, a galaxy may use up its gas by forming stars and thus have redder colors due 
to a lower star formation rate compared to the recent past. Alternatively, a galaxy's gas 
could be removed in denser environments and appear redder due to quenching of its star 
formation. Since we have quantified the environment of each galaxy (see §2.2), we are able 
to test whether this and other trends are due internal or external processes. 

The top right panel of Figure 4 suggests that both internal and external processes are 
important in controlling the trend between gas fraction and color. For gas-rich galaxies in our 
sample, with / gas > 0.4, the gas fractions are uncorrected with nearest neighbor distance. In 
contrast, galaxies with low gas fractions (/ gas < 0.4) are preferentially found within 0.5 Mpc 
of a luminous neighbor galaxy. Thus, it appears that some reduction in the gas fraction is 
related to internal processes such as normal star-formation, "primordial" removal of gas at 
reionization, or outflows, however, only environmental effects are very effective at removing 
all of the gas from a dwarf galaxy. This is further supported by the bottom right panel 
of Figure 4. While red galaxies only exist near luminous neighbors, the colors of the blue 
sequence galaxies do not depend on distance to the nearest neighbor. This result is similar 
to the relationship between color and environment, or morphology and environment, seen 
in more luminous galaxies (e.g. in Hubble 1936; Oemler 1974; Norberg et al. 2002; Blanton 
et al. 2005a), and the dwarf galaxy population of the Local Group (Grebel et al. 2003). We 
conclude that while internal processes are able to reduce the gas fractions to / gas ~ 0.4, 
external processes are require to lower the gas fractions (via removal of gas as suggested 
below via residuals in the Tully-Fisher relation) further. 



3.2. The Baryonic Tully-Fisher Relationship 

The Tully-Fisher relationship relates the mass of a galaxy as inferred from its luminous 
component to its dynamical mass, as measured by the maximum rotational velocity. The 
tight coupling between these components has greatly informed our understanding of galaxy 
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formation (Mo et al. 1998; van den Bosch 2000; Bullock et al. 2001; Mayer & Moore 2004). 
Dwarf galaxies with rotation velocities less than 90 km s -1 tend to fall below the optical 
Tully-Fisher relationship established by massive galaxies in the sense that the ratio of stellar 
to dynamical mass is lower for the dwarfs (McGaugh et al. 2000; Matthews et al. 1998). This 
"break" has sometimes been interpreted as evidence that dwarf galaxies have preferentially 
lost baryonic mass due to e.g. supernovae winds (Dekel & Silk 1986; Dekel & Woo 2003). 
McGaugh et al. (2000) first noted that this break is removed when both the stellar and 
gaseous components are taken into account. This baryonic Tully-Fisher relationship has been 
investigated by many groups (Bell & de Jong 2001; Verheijen 2001; Gurovich et al. 2004; 
McGaugh 2005). Understanding the baryonic Tully-Fisher relation for low mass galaxies 
should place tight constraints on the role of ram pressure stripping, tidal stripping, feedback, 
galactic winds, photoevaporation during reionization in the evolution of dwarf galaxies. 

In Figures 6 and 7, we explore the Tully-Fisher relationship within our dwarf galaxy sam- 
ple. We estimate the maximum rotational velocities for our galaxies based on the inclination- 
and turbulence-corrected HI profile half- width (W20^/2). While we will present resolved 
optical Ha rotation curves for this sample in a forthcoming paper, the HI line-widths are a 
more reliable measure of the maximum rotation velocity in these low mass systems. This 
is because the optical rotation curves of low mass galaxies rarely show a velocity turnover 
at the last optically measured data point (Matthews & Gallagher 2002), leading to an un- 
derestimate of the true maximum rotation velocity. Resolved HI rotation curves also do not 
always show a velocity turnover (McGaugh et al. 2001), however the HI tends to extend 
2-3 times beyond the observed optical radii (Begum et al. 2006) and thus is a more reliable 
probe of the true maximum rotation velocity. 

The dwarf galaxies in our sample have a significant rotation component (Figure 1). We 
correct the observed HI line-widths for line broadening due to turbulent velocity dispersion 
and inclination using the formula first proposed by Bottinelli et al. (1983): 

= W20-W20, 
sin % 

Where W20 is the observed HI line-width, W20 t is the turbulent velocity correction term 
and i is the inclination angle inferred from the optical images. We confirm the validity of a 
linear turbulence correction based on modeling the integrated velocity profiles of simulated 
galaxies constructed from Hernquist (1993) model disk galaxies. For nearby dwarf galaxies 
with rotation velocities similar to our sample, Begum et al. (2006) have measured velocity 
dispersion in the gas component of a\ os = 8 km s _1 based on the average line-of-sight velocity 
dispersion measured from 2D velocity maps. Using the Begum et al. value, this results in 
a turbulence correction of W20 t = 16 km s -1 , which we use here. Altering this turbulence 
correction (say to 25 km s~ x ) does not change our results below. 
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Using the approach of Haynes & Giovanelli (1984), we infer sini from the b/a axis ratio 
based on two-dimensional Sersic fitting of the r-band SDSS images. In particular, we assume 
the disks to be intrinsically oblate and axisymmetric, with a three-dimensional axis ratio of 
q , and take: 



with q = 0.19. Again, changing the particular value of q (say to 0.3) again does not change 
our results significantly. 

Figure 6 shows the stellar, HI and baryonic Tully-Fisher relationships for our dwarf 
galaxy sample. In the three panels of this figure, we linearly regress to fit a function of the 
form: 



where a and b are the fitted intercept and slope of the Tully-Fisher relationship, respectively. 
These fits are regressions with mass as the independent variable. Table 2 lists their values 
and uncertainties (calculated from bootstrap resampling). In addition, we calculate the 
correlation coefficient r between each pair of variables, listed in each panel of Figure 6. The 
correlation coefficient is zero for two unrelated quantities and unity for two perfectly linearly 
related quantities (e.g. Lupton 1993). 

The stellar Tully-Fisher (Figure 6, right panel) has significantly more scatter than the 
either the HI or total baryonic relationships (Figure 6, middle and left panels). The stellar 
relation has a much lower correlation coefficient (r ~ 0.3 versus r ~ 0.8) and considerably 
more scatter about the best fit regression line, as listed in Table 2 (about 50% scatter 
compared to about 30% scatter), than do the HI mass or baryonic mass TF relationships. 
This result is perhaps not surprising: at low luminosities, stars contribute far less to the 
total baryonic mass compared to massive galaxies (Figure 3). Thus, for dwarf galaxies of 
the same total mass, the amount of gas turned into stars varies between systems. While 
the scatter in the baryonic Tully-Fisher relationship is much smaller compared to the stellar 
relationship, it does have more scatter than accounted for by our observational error bars. 
The fitted slope is unchanged if we restrict our sample to only galaxies whose inclination 
corrections are small {b/a < 0.4) or if we alter our assumptions about turbulence or the 
intrinsic axis-ratio. 

In Figure 7, we compare our dwarf galaxy sample to Tully-Fisher data available in the 
literature. We compare only to literature samples for which HI and optical measurements 
are available: Haynes et al. (1999), Verheijen (2001), Matthews et al. (1998), and McGaugh 
et al. (2000). For the literature datasets, we recalculate distances based on H = 70Mpc 




(2) 




(3) 
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km s 1 and a model of the local velocity field ( Willick et al. 1997) (for galaxies with recession 
velocities less than 6000km s _1 ). For the Haynes et al. (1999) data, we assume a constant 
color ofV-I= 1.0. 

In the left panel of Figure 7 we plot the optical Tully-Fisher relationship in the J-band. 
McGaugh et al. (2000) first noted that low mass dwarf galaxies fall well below the optical 
Tully-Fisher (TF) relationship as defined by luminous spiral galaxies in the sense that, for 
a given rotation velocity, the dwarfs are optically under-luminous compared to the expected 
TF value. McGaugh et al. also noted that these same galaxies lie on the same relation 
defined by luminous galaxies in the baryonic TF relationship (Figure 7, right panel). Again, 
this is not surprising: we have shown above that dwarf galaxies have high gas fractions and 
the optical luminosity /stellar mass is often only a few percent of the total baryonic mass. 
However, unlike McGaugh et al., we do not see evidence for a distinct break in slope of the 
optical TF occurring at 90 km s _1 . At low luminosities, the scatter around the optical TF 
relation does increase considerably, but not preferentially to high or low velocities. 

That a single linear fit adequately describes the baryonic TF data over more than three 
decades of mass is significant. Many authors have suggested that supernovae-driven winds 
will preferentially remove baryonic material from galaxies with t> max < 100 km s^ 1 due to 
their far smaller gravitational potentials compared to luminous galaxies (Dekel & Silk 1986; 
Dekel & Woo 2003). Figure 7 suggests that a critical mass threshold does not exist below 
which mass-ejection is efficient down to baryonic masses of 10 8 M & and rotational velocities 
of as low as 20 km s _1 . If a critical threshold existed, we would have expected to see a break 
in the baryonic Tully-Fisher. 

While it might be possible to explain this single relationship via a baryonic mass-loss 
rate which varies smoothly with mass, the measured slope of this relation suggests that this 
is not the case. From the virial theorem, the expected slope of the baryonic TF is b = 0.33. 
The predicted slope from cosmological simulations is b ~ 0.29, or as more commonly quoted 
1/6 ~ 3.5 (Bullock et al. 2001; Kravtsov et al. 2004), in good agreement with our value of 
b = 0.27±0.01, or 1/b = 3.70±0.15. Thus, we interpret the baryonic TF relationship shown 
in Figure 7 as evidence that dwarf galaxies do not preferentially lose significant baryonic 
mass compared to luminous galaxies. 

Finally, we consider the residuals from the TF relationship. We define these residu- 
als in terms of the ratio W20 i:t /W20 ijt (TF), the corrected observed HI line- width divided 
by that expected from the fitted linear baryonic Tully-Fisher relationship. This quantity 
is defined in the sense that larger W20i tt /W20i tt {TF) values correspond to lower inferred 
baryon-to-dark matter ratios. Figure 8 shows this quantity as a function of projected near- 
est neighbor distance, color, gas fraction and axis ratio b/a. We find no correlation between 
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these residuals and the axis-ratio b/a, except that the very flat galaxies (b/a < 0.25) never 
exhibit low W20 ijt /W20 i:t (TF) ratios (Figure 8, bottom left panel). That no strong trend 
exists with axis-ratio suggests that the measured HI line-widths are properly sampling the 
galactic potential. Dwarf galaxies with low gas fractions have low W 20 j jt /W 20 i jt (TF) ra- 
tios. This suggests that low gas fraction dwarf galaxies, which are preferentially found in 
dense environments, have higher baryon-to-dark matter ratios compared to the majority of 
dwarfs in our sample. If, for example, reduced gas fractions in dense environments were due 
to interaction-induced star formation, we would expect these ratios to remain constant. If 
this were due to ram pressure stripping alone one might expect this ratio to be lower in 
dense environments: gas stripping would preferentially remove baryonic mass, reducing the 
ratio of baryons to dark matter. There is little correlation with color, indicating that if gas 
stripping is important, it is not directly effecting to the star-formation rates at the center. 
Alternatively, it is possible that this is simply a selection effect and the HI line-widths do not 
properly sample the galactic potentials of the lowest gas fraction galaxies. As the number 
of observed galaxies in this mass range increase, the TF residuals should place significant 
constraints on galaxy evolution processes important at these masses. 



4. Conclusions 

We present results from an HI survey of low luminosity dwarf galaxies in the absolute 
magnitude range —13.5 > M r — 51og 10 h 70 > —16 selected from the SDSS. These data rep- 
resent the largest homogeneous sample of low luminosity dwarf galaxies with well-measured 
HI and optical properties. The sample spans a range of total baryonic masses between 
Mbaxyonic = 10 8 - 10 9 M with measured HI line-widths between (v Iot ) = (W20 M /2) = 20-80 
km s _1 . The dwarf galaxies in this sample are found in a wide range of environments, from 
dense groups to truly isolated galaxies. 

For the majority of our dwarf galaxy sample, the total galactic baryonic mass is dom- 
inated by the gas mass rather than stellar mass. The average gas fraction for our sample 
is (/gas) =0.6, with gas fractions as high as 95%. This significantly exceeds that of gas-rich 
galaxies at higher luminosities. Thus, dwarf galaxies are inefficient at turning gas into stars 
compared to more luminous galaxies. We find that gas-rich galaxies have bluer g — r colors 
compared to gas-poor objects implying that gas-rich dwarf galaxies have had more recent 
star formation. The majority of galaxies in our sample are currently forming stars based on 
observed Ha emission in the SDSS spectra. We find that the current rate of star formation 
per unit area (based on the total w-band flux) is not correlated with gas fraction, confirming 
previous results suggesting that the star formation rates in dwarf galaxies are stochastic (van 



-16- 



Zee 2001; Hunter & Elmegreen 2004). Thus, dwarf galaxies form stars with a rate that varies 
over million year time scales, but is correlated with the gas fraction over longer time periods. 

Using the range of environments within our sample, we examine whether internal or 
external processes control the gas fractions of dwarf galaxies. Trends between galaxy prop- 
erties, such as color or morphology, and environment have been well studied for galaxies at 
brighter luminosities (Norberg et al. 2002; Blanton et al. 2005a), but have been previously 
limited to the Local Group and nearby groups at dwarf luminosities (Grebel et al. 2003). 
Our sample allows us to quantify these trends for dwarf galaxies over a much wider range of 
environments. We find that dwarf galaxies with gas fractions / gas < 0.4, comprising 15% of 
this sample, are exclusively found within 0.5 Mpc of a luminous neighbor galaxy. Galaxies 
with gas fractions / gas > 0.4 are found across the full range of environments contained within 
this sample. Thus, while internal galaxies processes can reduce the gas fractions, external 
processes are required to fully remove gas from a dwarf galaxy. 

The optical Tully-Fisher relationship within our dwarf galaxy sample demonstrates con- 
siderable scatter. As previous work has shown for brighter samples (Bell & de Jong 2001; 
Verheijen 2001; Gurovich et al. 2004; McGaugh 2005), the baryonic Tully-Fisher relationship 
within our sample is much tighter. When we consider the Tully-Fisher relationship across 
all absolute magnitudes, we find a logarithmic slope of b = 0.27±0.01 (or 1/b = 3.70±0.15). 
These results are in accord with CDM predictions of 1/b ~ 3.5 assuming no preferential 
loss of baryons as a function of mass (Bullock et al. 2001). Thus, our results suggest that 
processes which preferentially remove gas from dwarf galaxies are not important for the ma- 
jority of galaxies in our sample with maximum rotation velocities f max ~ 20 — 80 km s _1 
(roughly corresponding to total dynamical masses of M tota i = 10 9 — 1O 1O M ). 

If galaxies are associated one-to-one with ACDM dark matter halos, the mass-to-light 
ratios of galaxies must decline as a function of host halo mass (Benson et al. 2003). The most 
well-known example of this problem is the Local Group "substructure problem" (Klypin et al. 
1999; Moore et al. 1999). To explain this trend in mass-to-light ratio, investigators often 
invoke processes such as ram pressure stripping or supernova feedback that either removes gas 
or heats it sufficiently to prevent star-formation in low mass galaxies (e.g., Dekel & Silk 1986; 
Dekel & Woo 2003; Mori & Burkert 2000; Bullock et al. 2000). Most semi-analytic models 
require significant baryonic mass loss (Cole et al. 2000; Benson et al. 2003; Croton et al. 
2005), predicting baryon-to-dark matter ratios for dwarf galaxies well below the "Universal" 
baryon fraction (Read & Trentham 2005). Comparing the metallicity of galaxies over a 
wide range of baryonic masses, Tremonti et al. (2004) found that dwarfs are metal depleted 
relative to massive galaxies and argue in favor of galactic winds efficiently removing enriched 
material from low mass galaxies. However, at the masses of the galaxies in our sample, 
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detailed models of the physics of gas blow-out and blow-away tend to disfavor significant 
gas mass loss, predicting loss of only a few percent (Mac Low & Ferrara 1999; Ferrara & 
Tolstoy 2000). Our sample cannot rule out small amounts baryon mass loss from these 
galaxies, nor scenarios in which gas remains bound to the galaxy but is blown away from the 
region participating in star formation. This may explain the low observed metallicities of 
dwarf galaxies — enriched gas is removed from regions of active star formation, resulting in 
low observed metallicities, but either remains bound or represents only a few percent of the 
total baryonic mass. Indeed, recent cosmological simulations which include realistic feedback 
and chemical evolution suggest that the metallicity, mass-to-light ratios and other galaxies 
properties can be explained without including supernova energy feedback (Tassis et al. 2006). 
Our results favor such dwarf galaxy formation models in which processes such as supernovae 
winds are not responsible for significant gas or baryon mass in low mass galaxies. 
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Table 1. Low Mass SDSS Dwarf Galaxy Properties 
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Table 1 — Continued 



Name 


a (J2000) 


5 (J2000) 


Distance 


m r 


M r 


9 - r 


Moff ,r 


r ott,r 


b/a 


M stellar 


m H i 




W20/2 


W20 it /2 


CT W20 i>t /2 




(h : m : s) 


(":':") 


(h~Q Mpc) 










kpc (") 




(logMQ) 


(logMQ) 


(LogMQ) 


km s~ 1 


km s~ 1 


km s - 1 


421787 


23 31 02.1 


-00 07 55.0 


24.1 


17.16 


-15.63 


0.36 


21.7 


1.0 (8.5) 


0.29 


8.02 


8.60 


0.05 


61.6 


54.9 


2.3 



Note. — Galaxy names are based on position in the Blanton ct al. (2005) low luminosity SDSS catalog. Absolute magnitudes and g — r colors have been corrected for Galactic extinction based 
on Schlegel et al. (1998). Distances are estimated based on a model of the local velocity field (Willick et al. 1997) and assume Hq = 70hjQ Mpc km s — 1 . The effective surface brightness (Meff,r)i 
effective half-light radius (r c ff r ) and axis ratio (b/a) are measured from the SDSS r-band images. Stellar masses are calculated based on Bell et al. (2003) and a Kroupa ct al. (1993) initial mass 
function; HI masses arc based on our Arecibo and GBT observations. The HI velocity half-width, W20/2, is the observed uncorrected 20% peak HI half-width. The quantity W20j f/2 is the HI 
half-width corrected for inclination and turbulence as described in §3.2. Errors on the HI mass (cjvf HI ) and HI line-width (<^\Y20- f/2) wcrc calculated from the scatter in the mean quantities 
recovered from Monte Carlo simulations. Galaxies which wcrc not detected in the radio do not have HI properties listed in this table. 
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Table 2. Tully-Fisher Relation Fit Parameters 



Galaxy Sample 


a 




/ 






a 


Stellar mass (current sample) 


-0.14 


± 


0.02 


0.43 


± 


0.06 


0.12 


Baryonic mass (current sample) 


-0.32 


± 


0.05 


0.53 


± 


0.08 


0.13 


Mj (full sample) 


0.37 


± 


0.01 


-0.09 


± 


0.01 


0.11 


Baryonic mass (full sample) 


-0.19 


± 


0.02 


0.27 


± 


0.01 


0.10 



Note. — Tully-Fisher linear regression fit parameters for Equation 3. 
The zero-point a is in units log 10 (W20j i t/2/(50 km s -1 )). The slope b is 
in units log 10 (mass/10 8 /if 2 M o ) per log 10 (W20 M /2/(50 km s" 1 )). In the 
J-band Tully-Fisher relationship we replace log lo (mass/lO 8 /iy O 2 M ) with 
Mi — 51og 10 /i7o + 20. The quantity a represents the logarithmic (base 
10) scatter around each fit. The first two rows ('current sample') are 
fits to the galaxy data presented in this paper; the last two rows ('full 
sample') are fits to this sample plus data from the literature as desribed 
in §3.2. 
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Fig. 1. — SDSS r-band images and HI line- widths for three representative dwarf galaxies. 
The images are V on a side. We list the absolute r-band magnitude and optically measured 
axis ratio (b/a) for each galaxy. Of the dwarf galaxy sample presented here, 18% have double- 
peaked HI profiles similar to the top panel, 26% have flat-topped HI profiles, and 56% are 
single peaked profiles. If we consider only edge-on galaxies (b/a < 0.4), the percentage of 
double-peaked profiles increases to 30%. The distribution of HI profile shapes suggests that 
the majority of dwarf galaxies in our sample have regular rotation fields. 
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Fig. 2. — The distribution (top left to bottom right) of absolute magnitude, distance, effective 
surface brightness, effective half-light radius, Sersic index and axis ratio for our dwarf galaxy 
sample. All photometric parameters are measured from the SDSS r-band images. Distances 
are estimated based on a model of the local velocity field (Willick et al. 1997). In each panel, 
we list the average value of each quantity for the 101 dwarf galaxies in this sample. 
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Fig. 3. — The gas fraction, / gas = M gas / (M gas + M ste n ar ), plotted as a function of absolute r- 
band magnitude. The majority of dwarf galaxies in this paper (black circles) are dominated 
by the gas mass, rather than stellar mass, with gas fractions as high as 95%. Galaxies 
which were not detected in HI are plotted as upper limits. For comparison, we plot data 
from Schombert et al. (2001, squares), Bell & de Jong (2000, upward triangles) and Haynes 
et al. (1999, downward triangles). We infer from this plot that, compared to more luminous 
galaxies, dwarf galaxies are far less efficient at turning gas into stars over their lifetimes. 
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Fig. 4. — Top Left: The gas fraction, / gas , plotted as a function of g — r color. Assuming 
the g — r color traces recent star-formation activity, the fraction of star-formation occur- 
ring in the recent past is higher for the gas-rich galaxies in our sample. Bottom Left: In 
contrast, the current star formation rate per unit area, as measured by the SDSS -u-band 
flux, is not correlated with gas fraction. Top Right: The gas fraction plotted against the 
projected distance to the nearest luminous neighbor. Galaxies with gas fractions / gas < 0.4 
are preferentially found within 0.5 Mpc of a luminous neighbor galaxy. The gas fractions of 
galaxies with / gas > 0.4 are uncorrelated with nearest neighbor distance. Bottom Right: We 
also compare the g — r color and nearest neighbor distance distribution of our sample to that 
of the parent SDSS dwarf galaxy catalog {small black symbols). 
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Fig. 5. — Left: The stellar birthrate timescale defined as M ste iiar/<S'F-R to t. The median 
birthrate timescale for this dwarf sample is 15Gyr (black arrow). Assuming a single galaxy 
formation time of ~12Gyr (grey arrow), roughly half the galaxies in our sample had lower 
star formation rates and half have had higher rates in the past than presently observed. This 
suggests that star formation in dwarfs is stochastic and that the current star formation rate 
is not necessarily representative of the time-average rate. Right: The large gas consumption 
timescale (M gas / SFRt ot ) suggests that the majority of galaxies will be able to make stars at 
the current rates for another Hubble time. 
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Fig. 6. — The stellar mass, HI mass, and total baryonic mass plotted as a function of 
rotation velocity as measured by the inclination- and turbulence-corrected HI half line- widths 
(W20j )t /2). The solid line is the (inverse) regression fit described in §3.2, whose parameters 
are quoted in Table 2. Each panel lists the correlation coefficient r between the two quantities; 
the correlation coefficient is zero for two unrelated quantities and unity for two perfectly 
linearly related quantities. The scatter in the stellar mass relationship (left) is much larger 
compared to the HI (middle) and baryonic Tully-Fisher relationship (right). 
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Fig. 7. — The optical and baryonic Tully-Fisher relationships: absolute /-band luminosity 
(left) and the total baryonic mass (right) plotted as a function of maximum rotational velocity 
(W20i ) t/2). We compare our dwarf galaxy sample to literature data for which HI and optical 
measurements are available: Haynes et al. (1999), Verheijen (2001), Matthews et al. (1998), 
and McGaugh et al. (2000). The line in each panel is the fit we describe in the text. The 
slope in the right-hand panel b = 0.27 ±0.01 is consistent with the CDM prediction b ~ 0.29 
(Bullock et al. 2001). 
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Fig. 8. — Residuals from the baryonic Tully-Fisher relationship, defined as the ratio 
W20 itt /W20 itt (TF) of the actual line width to the expected line width from the Tully-Fisher 
relationship. We show this ratio versus nearest neighbor distance, g — r color, gas fraction 
/g as and optical axis ratio b/a. 



